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Abstract 

The frequency and shear dependent critical viscosity at a correla- 
tion length ^ = K,~^, has the form r] = r]QK~^'^G{zi, Z2), where zi and 
Z2 are the independent dimensionless numbers in the problem defined 
as zi = aiid Z2 = ■ "^^^ decay rate of critical fluctuations of 

correlation length is Tqk^ and is the effective wave number for 
which Tok^ = S, the shear rate. The function G{zi,Z2) is calculated 
in a one loop self-consistent theory. 

PACS number (s):64.60Ht 
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Response functions in general tend to diverge near a second order phase 
transition point. As the fluctuations become more long ranged and drive the 
system towards the transition, the response functions in the thermodynamic 
limit (if static) and in the hydrodynamic limit (if dynamic) diverge. Some 
like the susceptibility (static) and thermal diffusivity (dynamic) [1] diverge 
strongly with an exponent one, others like the specific heat (static) and 
shear viscosity (dynamic) [2-4] diverge weakly with an exponent close to 
zero. The small exponents have always provided an interesting challenge 
for theorists and experimentalists. The shear viscosity of ordinary fluids or 
binary mixtures, in particular, has been a favourite candidate for pushing 
theory and experiment to the limit. The exponent is small but over the 
last two decades the experiments [5,6] have increased the accuracy of the 
measurements so strongly that theorists have had to worry about higher 
loops [7,8], something which is rarely done. 

The critical divergence is masked if the system is not in the hydrodynamic 
limit. Hydrodynamic limit implies the wavelength is larger than the critical 
correlation length and all frequencies are lower than the rate of decay of 
critical fluctuations. While wavelengths are always larger than the critical 
correlation length, critical slowing down may imply that frequencies may 
not be smaller than the fluctuation relaxation rate. The frequency of an 
oscillating viscometer used to measure the critical viscosity can in that case 
determine the viscosity of the sample and the fluid will become viscoelastic. 

The frequency dependence of the viscosity predicted [9,10] and measured 
over the last two decades. It has been recognized by Oxtoby that a shear 
rate will also introduce a time scale the problem and if the fluctuations are 
more long lived than this time scale, then experiments on this time scale 
will not feel the longest-lived fluctuations and the viscosity will be limited 
by the shear rate - a phenomenon known shear thinning. The dependence 
of various static quantities on the shear rate as well as its effect on light 
scattering has been considered by Onuki [12] as well as Onuki and Kawasaki 
[13]. With so much information available on both theoretical and experi- 
mental fronts, experiments have been designed to test the shear thinning of 
critical viscoelastic fluid. We have carried out a one loop calculation of the 
frequency and shear dependent critical viscosity and present our results in 
this communication. 

In the hydrodynamic regime, characterized by the sole length scale - the 
correlation length ^ = = ^o[^^^^]~'') with i/ ~ |, the shear viscosity 
diverges as 

r]{K) = TyoK"^" (1) 

where x,; is the small exponent discussed above which is found to be around 
0.068. The decay time of the fluctuations, r , is given by 

T-' = r(k,K)^-k^ (2) 

X 

where L is the Onsager coefficient and % is the static susceptibility. The 
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critical viscosity can , to a very good accuracy , be taken as x^^ik,^) = 
k'^ + and the Onsager coefficient or thermal diffusivity diverges at the 
critical point. The behavior of wavelength and correlation length dependent 
diffusivity is governed by the Kawasaki function [1] which was simplified for 
practical use by Ferrell and can be used in the form 



L = ro(r + k')-^ (3) 

in D = 3, where D is the spatial dimensionality. The characteristic decay 
time of a fiuctuation of wavelength equal to the correlation length is given 
by 

r-^ = T{k) = ToK^ (4) 

Hydrodynamic regime implies that frequencies are such that cur -C 1. How- 
ever, at a fixed frequency u>, as k. decreases on approaching the critical point, 
r diverges and it is not possible to satisfy lot < 1. In that situation Eq.(l) 
can not hold and if ~ 0, the response is limited by the frequency u. 
Since u) scales as k^, it is clear that the limiting viscosity will be of the form 

r]{uj) = rjo{-iuj)~^ (5) 

and the full viscosity will be governed by 

r]{K,(jj) ^r]oK-''^F{z)-^ (6) 

where F(z) is a function of the dimensionless variable z = t^^^- For z — >■ 
0, F{z) — > 1 and if 2; 3> 1, F{z) oc z. The simplest possible functional 
representation of F{z) is 

F{z) ^1 + Pz (7) 

Now, /3 is a number of 0(1), which can be found from the small z form of 
the one loop integral or from the z ^ 1 form. If determined from the low 
frequency end /9 = || ~ 0.59. From the high frequency end /? ~ 0.2. The 
two loop calculation is found to enhance /? by about 30% to 0.8 at the low 
frequency end. The experimental value found by Berg et al is about 1.2. 

If a shear rate S is now switched on which result in a mean fiow in (say) 
the x-direction, then the velocity can be written as 

v = Sye^ (8) 

The shear rate S introduces a new length scale defined by 

S = Tkl (9) 

Strong shear implies kc > k, while the reverse is the case of weak shear. 
Our primary interest will be in strong shear which will always be the case, 
sufficiently close to the critical point. For the shear thinning of a viscoelastic 
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fluid we now have two dimensionless numbers in the problem Zi — and 
Z2 = 2rpfc3 s-iid Eq.(6) is generalized to 

r)^rjoK-^-[G{z,,Z2)]-'^ (10) 
The one loop answer for G{zi,Z2) is 

G^{z,,Z2)^[l + 0.6z^][l + ^x--4-? (11) 

2^1 ^ + - 
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in the case of strong shear i.e. kc > k, which in the language of Zi and Z2 
means ^ > 1. This is the one loop result. 

Using the fact that the one loop scale factor viscoelasticity is small by a 
factor of two compared to the experiment, one can use the improved formula 



G(^i, Z2) = [1 + 1.2^i][l + A X j^]' (12) 

2^2 zi 

The central results of our work are contained in Eqs.(lO) and (12). We now 
sketch the derivation . The equation of motion for the order parameter is 
nonlinear. The effect of the nonlinearity is to make the transport coefficient 
divergent in the absence of the shear. The contribution of the shear to the 
equation of motion is a linear term. Consequently, we will work with an 
effective equation of motion where the effect of the nonlinear terms will be 
handled by a dressing of the transport coefficient. Consequently in momen- 
tum space, the order parameter (/){x) satisfies 

where x"^ — k"^ + and L — Tqx^ as explained in Eq.(3). We now need to 
work out the susceptibility in the presence of S and this is simply retracing 
Onuki's calculation with the present L and this leads in a straightforward 
fashion to 

- k' + + kii^-^-y^ (14) 

The effect of the changed diffusion coefficient shows up in the slightly different 
kx dependence in the last term on the right hand side. In the expression for 
viscosity, the susceptibility will be averaged over all directions and for that 
specific case, we will use an angle averaged form of x^^, where we replace 
\kx\'^ by A; 2 < sin^Ocos^O >= |A;2 . The one loop shear viscosity is now given 

by 

■n{K,uj,S) = — / d p Y 

47ry '^-iu; + 2ToP^[x{p,K,S)]-2 
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''[^i+P^b^ + i(f)^(t)M + i]^] 



(15) 



It is the characterization of this one loop integral which has to be carried 
out. The limit Z2 ^ oo (or 0, i.e. no shear) has already been treated 

and this is what we represent to the lowest order by Eq.(7). It should be 
noted that Eq.(15) has a logarithmic divergence rather than a power law 
divergence. To extract G{zi,Z2) from Eq.(15), we use the fact that ^ 1 
and expand Eq.(lO) as 

rj = r]oH-^'^[G{zi,Z2)]~^ 

X 

~ ?7o[l - Xj^lriK - -^lnG{zi, Z2)] 
^ r)B + Voiin— - ]-lnG{zi, Z2)] 

K 6 

(16) 

where tjb is a background viscosity and the part within the square brackets 
emerges from the loop calculation. 
For Z2 — > 00, Eq.(15) becomes 

r]{uj, K) = —— / dp 



2roio '^{l+pY[z,+p^{l+p-^)l] 
nn^ + ln2 - 2+ ^-^z, + ....] 



2ro' K 4 16 



(17) 



for zi -C 1. On the other hand for Zi ^ 1, 77(0;, n) — ^[In^ — \lnzi + ], 

which when combined with Eq.(17), leads to the approximation shown in 
Eq.(7). 

Another limit which can be similarly explored is a; — > 0, i.e. 2:1 —> and 
Z2 — > 0, but ^ 7^ in which case, Eq.(15) becomes 

riiS, = ^ dp f (18) 

For ^ <g 1, an expansion similar to Eq.(17) obtain. For ^ » 1, we scale 
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momenta by (— )3 to find 



2roio [p2 + + 



2ro^ kc 3 



Combining Eq.(19) with the small ^ form, 



?7o „ A 4, ^ 8 27r 



Gizi ^ 0, ^2 ^ 0) = [1 + 0.3(-)5]2 (20) 
Finally as re — > 0, 2:1 — > oo with finite and Eq.(15) reduces to 

s)^^ dp , ^ , I 1.,. — ^ , — r^-TT (21) 

For Z2 — > 0, we recover Eq.(19), while for Z2 ^ 00 

^(^, S) = ^[inA - l/nz2 - ^ + ...] (22) 

We thus obtain the following limiting forms: 

G{zi, Z2) — > 1 + O.Gzi i/ Z2 — > 00 

3 1 

G{zi,Z2) 0.62;i[l + -— ]^ if zi^ 00 

z ^2 



^1 

G'(^i,^2) ^ [l + 0.3(^)^]2 ^1,^2^0 

^2 



(19) 



(23) 



The form of G{zi, Z2) taking into account all the above constraints is given by 
Eq.(ll) and with the phenomenological improvement, we arrive at Eq.(12). 
Very recently Berg [15] has arrived at a phenomenological rule for shear 
thinning. The principle difference between Eq.(ll) and Bergs' observations 
lie in the existence of the nonanalytic behavior in G{zi,Z2)- This could be 
an interesting issue during comparison with experimental data. 
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